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c/3 \ Making use of the tight-binding Green's function technique, we have calcu- 

lated the STM images of graphite with surface and sub-surface defects, while 
c5 \ taking into account the relaxation of the lattice due to defects. We have 

__j ■ demonstrated that two different physical mechanisms may result in the for- 

a, 

O \ mation of hillocks in the STM images: buckling of the graphite surface due to 
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interstitials between the uppermost graphite layers and the enhancement of 
the electron density of states close to the Fermi energy on the carbon atoms 
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O ■ in the vicinity of vacancies. Our results indicate that small hillocks may orig- 

0> . inate both from the interstitial clusters and from the vacancies. By contrast. 
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H . defects 



however, large hillocks in excess of 10 A in diameter can be caused only by 
interstitial clusters. 
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Scanning tunneling microscopy (STM) is a powerful tool for investigating surfaces with 
defects inasmuch as it gives real-space images of surfaces with atomic resolution. Radia- 
tional defects near the surface profoundly effect STM images of graphite. Appearances of 
hillock-like features in STM images of graphite after irradiation with noble gas ions have 
been reported repeatedly []I], 0, ^ ^ ^, ^. 

These hillocks (or bumps or protrusions), which reveal themselves as a substantial 
increase in the tunneling current, have an average lateral diameter of from several A up 
to tens of A, with the typical height of hillocks constituting 1-5 A. Since these structures 
have never been observed on graphite surfaces free of defects, they are associated with 
defects created by incident ions near the surface. 

Despite numerous studies, however, the actual nature of hillock-like features remains 
controversial. Hillocks may result from local buckling of the surface due to internal stresses 
induced by subsurface defects ||l| and- -more specifically- -by carbon interstitial clusters 
between the graphitic layers |0, 0, i.e. within the framework of this assumption the 
hillocks can be interpreted topographically. 

Recent molecular dynamics (MD) simulations 0, § provide evidence in favor of this 
interpretation. It has been demonstrated that ion irradiation of graphite may result in 
the formation of interstitial clusters under the graphite surface and in the buckling of the 
surface layer. 

The physical mechanism for interstitial cluster formation is proposed in |^ . The inter- 
stitials created by ion impact lie predominantly between graphite layers and have a 
high mobility parallel to the layers. At the same time, as has been shown in [§] by the 
MD method, an interstitial stretches the graphite lattice within its vicinity, with other 
interstitials in that same interlayer region tending to be attracted to the stretched region. 
Thus, due to the minimization of the total buckling of the graphite layers, the formation of 



interstitial clusters lowers the deformation energy, as compared to the case of interstitials 
being spatially separated. 

A different interpretation [Q of hillock-shaped features is that the surface atomic va- 
cancy is seen as a protrusion in STM images of graphite because of the vacancy-induced 
enhancement of an electron density corresponding to the states close to the Fermi en- 
ergy (Ep) on the carbon atoms near the vacancy. Inasmuch as STM probes only some 
of the electronic states of the sample within a certain energy range near Ep determined 
by the bias voltage Vuas, the enhancement of electron density leads to an increase in the 
tunneling current. Real geometry of the surface remains practically flat. 

It should also be noted that the hillocks may be explained by incident ions trapped 
near the graphite surface |0, but, in many cases (e.g., in the case of incident ions with 
high energy), the formation of hillocks seems to be associated with carbon atoms only. 
Adsorbate effects are believed Q to be either absent or relatively negligible. Therefore, 
we restrict our consideration only to cases of carbon interstitials and vacancies. 

The lack of a clear understanding of the nature of the hillocks may be due in part to 
the comparatively small number of theoretical works on the subject. The effect of lattice 
vacancies on the STM image of graphite was studied in |^ by the tight-binding (TB) 
molecular orbital method and in [|1^ using the first-principles plane wave formalism. The 



results of simulations support the interpretation that vacancies may lead to the formation 
of hillocks. 

However, we are not aware of any theoretical works dealing with the simulation of the 
STM images of graphite in the presence of interstitials near the surface. 

In this paper, we simulate the STM images of the graphite surface having defects such 
as interstitials and vacancies making use of Green's function calculations based on the TB 
model. Our main goal is to clarify whether the vacancies and interstitial clusters under 



the graphite surface result in the formation of protrusions in the STM images. 

In general, in our simulations we adopt the technique used in |TT]. To describe the 



graphite electronic structure, we employ the parametrization of the TB Hamiltonian for 



carbon systems proposed in [|12[. Since the density of states near Ep in graphite is 
determined by 7r-states perpendicular to the graphite layers [|n|, |T^, we take into con- 
sideration only 2pz electrons. The tip was modeled as the final atom of a semi-infinite, 
one-dimensional chain. To the first order in the tip-surface interaction, the tunneling 
current J as a function of the tip coordinates in the (x, y) plane (parallel to the graphite 



surface) and the tip height h may be written as follows [11 
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where V is the tip-surface matrix element, the sum over i runs over graphite surface sites 
effected by the tip, Psurf{E) and ptip{E) are the local densities of states (LDOS) of the 
non- interacting surface and tip, respectively. 

To account for the deformation of the graphite lattice near the defects, the MD method 
was employed. Details of calculations are given in 0. Graphite was modeled by the 
crystallite consisting of six graphite layers, each composed of 1096 atoms. Crystallites 
with different numbers of atoms were also used, but the results proved to be practically 
independent of the crystallite size. Interstitials were initially distributed randomly over 
the crystallite between the uppermost layers. Ten initial variants were considered for each 
number of interstitials. The tunneling current was calculated after the relaxation of the 
lattice and interstitial cluster formation due to deformation interaction. 

We have calculated the STM images of the graphite surface near both a single interstitial 
as well as clusters composed of 2-10 interstitials. Fig.l (a) shows an isometric plot of the 
variation in the tip height h for a scan across (x, y) near ten- interstitial cluster on the 
graphite surface at a current of 4 nA and Vbtas = 0.08 V. Fig. 1(b) depicts a contour map 



corresponding to Fig. 1(a) and simulating the STM image. A dramatic protrusion above 
the interstitial cluster is evident. This hillock is due to a graphite surface buckling(?^ 1 A) 
which results in a corresponding growth of V and in an increase in the tunneling current. 
The graphite LDOS does not change substantially near the cluster. The positions of 
maxima and minima of the tip height near the cluster correspond to those in the defect- 
free case (in graphite every second atom is observable due to electronic effects [|13|)- Slight 
anisotropy in the corrugation (^ 0.05 A ) near the cluster is associated with the non- 
uniform distribution of electron density because of the anisotropic form of the interstitial 
cluster, see Fig. 1(b). 

Hillock-like features were observed for all clusters, but the height and lateral diameter 
of protrusions were dependent on the number of interstitials in the cluster Ni^t as well 
as on cluster shape. The average diameter of hillocks d as a function of Nint is shown in 
Fig. 2. It is seen that d increases linearly with Ni^t- The average height also grows with 
Nint, but rather slowly (h ^ 1 A for Nint = 1 and h ^ 1.5 A for Nint = 10), which may 
be associated with the two-dimensional character of interstitial clusters. 

Let us proceed to the effect of vacancies on the STM image. Fig. 3 shows the variation 
in the tip height h for a scan on the graphite surface near a single vacancy positioned on 
the (3 surface site (there is no atom below the /3 site in the second graphite layer). The 
tunneling current and Vbias were the same as in the case of interstitials. As follows from 
Fig. 3 and although real geometry remains flat, the vacancy results in the formation of 
a hillock with the height of ~ 0.8 A. The increase in h is governed by the enhancement 
of the LDOS at Ep on the atoms in the immediate vicinity of the vacancy, which is in 
agreement with the results of |^. 

The superstructure (SS) ^/S x ^/3i?30° is seen in the image. The SS was also evident 
in our calculated images near single interstitials. Similar SS were experimentally imaged 



next to some hillocks in |l|]. Our calculations confirm the interpretation [|I4| that these SS 
are actually a purely electronic effect due to point defects and are not indicative of any 
surface reconstruction. Note that SS are absent in Fig.l since interstitial clusters result 
in a long-range perturbation of the lattice. 

Qualitatively similar results were obtained for the vacancies in the a positions, diva- 
cancies, and trivacancies. The lateral size of the vacancy-induced hillock depends slightly 
on the number of vacancies in the system, see Fig. 2. Therefore, our results indicate that 
large, experimentally-observed hillocks exceeding 10 A in diameter can be caused only 
by the interstitial clusters. 

In summary, we have calculated the STM images of graphite having surface and sub- 
surface defects. Two different physical mechanisms may give rise to the appearance of 
hillocks in the STM images: buckling of the graphite surface due to interstitials between 
the uppermost graphite sheets and the enhancement of the partial charge density of states 
on the carbon atoms near vacancies. Small hillocks (< 10 A in diameter) may originate 
from both the interstitial clusters and vacancies. However, large hillocks (> 10 A) may 
be caused only by interstitial clusters. Thus, our results confirm the mechanism [§ of 
interstitial cluster formation in the irradiated graphite due to deformation interaction. 
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Figure captions 

Fig.l. Constant current mode of operation of STM: (a) Variation in tip lieiglit h for a 
scan across x, y on tlie grapliite surface near a ten-interstitial cluster at a current of 4nA 
and Vbias = 0.08 V. (b) The contour map of 1(a). Dots correspond to positions of carbon 
atoms in the surface layer, squares correspond to carbon atoms in the interstitial cluster. 

Fig. 2. The average diameter d of the hillock as a function of the number of interstitials 
and vacancies in the system. 

Fig. 3. The same as in Fig.l, but for a single vacancy. 



This figure "figl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/9910498vl 
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This figure "fig3.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/9910498vl 



